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ABSTRACT 

The growth history of large-scale structure in the Universe is a powerful probe of the 
cosmological model, including the nature of dark energy. We study the growth rate 
of cosmic structure to redshift z = 0.9 using more than 162,000 galaxy redshifts from 
the WiggleZ Dark Energy Survey. We divide the data into four redshift slices with 
effective redshifts z = [0.2, 0.4, 0.6, 0.76] and in each of the samples measure and model 
the 2-point galaxy correlation function in parallel and transverse directions to the line- 
of-sight. After simultaneously fitting for the galaxy bias factor we recover values for the 
cosmic growth rate which are consistent with our assumed ACDM input cosmological 
model, with an accuracy of around 20% in each redshift slice. We investigate the 
sensitivity of our results to the details of the assumed model and the range of physical 
scales fitted, making close comparison with a set of N-body simulations for calibration. 
Our measurements are consistent with an independent power-spectrum analysis of a 
similar dataset, demonstrating that the results are not driven by systematic errors. We 
determine the pairwise velocity dispersion of the sample in a non-parametric manner, 
showing that it systematically increases with decreasing redshift, and investigate the 
Alcock-Paczynski effects of changing the assumed fiducial model on the results. Our 
techniques should prove useful for current and future galaxy surveys mapping the 
growth rate of structure using the 2-dimensional correlation function. 
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1 INTRODUCTION 

The clustering pattern of galaxies, which is driven by the 
competing forces of gravitational attraction and universal 



cosmic expansion, is one of the most important probes of 
the cosmological model. In particular, measurements of the 
growth rate of cosmic structure, inferred from the patterns 
imprinted in the galaxy clustering distribution, can poten- 
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tially distinguish between the increasing number of alterna- 
tive theories of gravity which currently compete with Gen- 
eral Relativity to explain the observed acceleration of the 
universal expansion without the need for an exotic compo- 
nent like Dark Energy (see Tsujikawa 2010 for a review). In 



general the cosmic growth rate will evolve differently in these 
models, even if they have the same expansion history, allow- 
ing constraints to be placed on different theories of gravity 



(Linder & Jenkinsl 20031 iLinder & Cahn 2007 Guzzo et al 



2008 Wang 2008 



It has long been known in standard gravity that the 
growth rate / = d In D/d In a (in terms of the linear growth 
factor D and cosmic scale factor a) is well- approximated by 
f(z) ~ Q m (z) 0-6 , where £l m (z) is the matter density rel- 
ative to critical density at redshift z (Peebles| |1980|). More 



recent theoretical investigations ( [Linder &; Cahn |2007| have 



generalized this equation to a phenomenological relation 

f( z ) = n m (zy (i) 

where the parameter 7, the "gravitational growth index" 
will take on different values depending on the considered 
theory of gravity, where 7 ~ 0.55 for General Relativity in 
a Universe dominated by a cosmological constant A. In this 
study we test the self-consistency of the ACDM cosmolog- 
ical model on large scales, by measuring the growth rate 
from a galaxy survey assuming a fiducial ACDM model and 
comparing these measurements with the predicted values. 

A powerful technique for determining the growth rate is 
to measure and model redshift-space distortions imprinted 
in the 2-point galaxy correlation function (hereafter £ s (cr, tt)) 
of a galaxy redshift survey, where the two dimensions (a, tt) 
correspond to galaxy pair separations transverse and par- 
allel to the line-of-sight. This is possible because the bulk 
flows of matter produced by large-scale structure formation 
leave a pattern of correlated peculiar velocities in the red- 
shifts measured by galaxy surveys, introducing a character- 
istic distortion into the shape of the measured £ s (cr, tt). 

The models used to fit the distorted pattern present in 
£ s (cr, 7r) have their origins in the works of Kaiser (1987) and 
Hamilton (1992}. The principal ingredients which make up 



these models are: 

• The underlying real-space correlation function (£ r ), 
which describes the clustering of galaxies in the absence of 
redshift-space distortions. This function cannot be measured 
directly from real observations due to galaxy peculiar veloc- 
ities, but can be modelled from theory or simulations. 

• A prescription for the large-scale linear distortion of the 
redshift-space correlation function due to galaxy bulk flows. 

• A function for the galaxy pairwise velocity distribution, 
F(v), which represents non- linear effects in the velocity field 
resulting from small-scale interactions between galaxies and 
their surroundings. An example of these effects is the char- 
acteristic spike known as the "finger-of-god" observed in 
£ s (cr, 7r) at small tangential separations (a < 4/i _1 Mpc), 
due to the virialized motions of galaxies in dark matter 
haloes. 

• A prescription for the bias with which galaxies trace 
the underlying matter density field. 

Various approaches have been used in the literature for 
modelling these different components. In the simplest ver- 
sion of the model a power law is assumed for £ r on small 



scales, an exponential function for F(v), and galaxies are 
assumed to trace the matter-density field with linear bias on 
large scales. This simple model has been applied to analyses 
of many galaxy surveys including recent large spectroscopic 
surveys such as the 2-degree Field Galaxy Redshift Survey 



(2dFGRS, Colless et al. 2001) and the Sloan Digital Sky 



Survey (SDSS, York et al. 2000) 



The 2dFGRS spans an area « 1500 deg mainly located 



in two regions of sky. I Hawkins et al. (2003) used 166,000 



2dFGRS galaxies in the redshift range z — [0.0, 0.2] to mea- 
sure £ s (cr, tt) up to separations of 30/i _1 Mpc, constraining 
the growth rate with an accuracy of approximately 10% at 
an effective redshift z = 0.15. The SDSS Luminous Red 
Galaxy (LRG) sample contains more than 75,000 objects in 
the more extended redshift range z = [0.15, 0.47] spanning a 
volume bigger than 1 h~ 3 Gpc 3 . Its correlation function al- 
lows a measurement of / at an effective redshift of z — 0.35 



flCabre fe Gaztanaga[ [2009] |Okumura et al.||2008||Samushia 



et al. 2012). Noisier measurements of the growth rate have 



also been determined from higher-redshift surveys such as 
the VVDS (V IMOS VLT Deep Survey) at redshift « 0.8 
(IGuzzo et al. 1 120081). 



Other analyses of redshift-space distortions in the 
galaxy distribution, involving a wider range of cosmologi- 
cal parameters, include: the studies of Percival et al. (2004) 



who model the redshift-space distortions of the 2dFGRS by 
expanding its density field in spherical harmonics; Tegmark 
et al. (2006) who measure galaxy power spectra from SDSS 
LRGs at z ~ 0.3 and fit for a large set of cosmological param- 
eters; and da Angela et al.| ( |2008| who use QSO catalogues 
from the 2QZ and 2SLAQ surveys to compute and model 
£ s (cr, tt) for z <^> 1.5. These studies generally found agree- 
ment between the measurements and the predictions of the 
prevailing ACDM cosmological model. 

Although galaxy redshift surveys currently provide the 
most accurate means of measuring the cosmic growth his- 
tory, other methods are also being actively developed to ex- 
tract this information. These include studies of the luminos- 
ity function of X-ray clusters and their physical properties 
such as gas mass fraction ( |Rapetti et al.[ |2009| ), measure- 



ments of bulk flows in the local neighbourhood (Abate & 



Lahav 



2008 



Watkins et al. 



2009 



Nusser & Davis 



2011), 



and weak lensing distortions mapped in large, high-quality 



photometric imaging surveys (Heavens et al. 
2009] ). 



stein, 



2007 Bern 



In this study we analyze redshift-space distortions in 
the 2-point correlation function of the WiggleZ Dark En- 



ergy Survey (hereafter WiggleZ, |Drinkwater et aL] (2010)). 
This is a new high-redshift survey which has observed more 
than 200,000 galaxy redshifts over ~ 1000 deg 2 of sky in the 
redshift range z < 1. The availability of the WiggleZ data 
enables accurate cosmological measurements to be carried 
out for the first time in the high-redshift (z > 0.4) Uni- 
verse, the epoch at which the Universe apparently started 
its accelerating expansion. WiggleZ galaxies are selected as 



star-forming (blue) galaxies by photometric criteria (Blake 
et al. 2010|[Drinkwater et al. 2010| ), in contrast to previous 



surveys wich targetted Luminous Red Galaxies, and there- 
fore also allow us to investigate if the cosmological conclu- 
sions depend on the type of galaxy used to trace large-scale 
structure. 

The growth rate has been measured from the WiggleZ 
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dataset by Blake et al. (2011a) using a power-spectrum anal- 



ysis. In this study we carry out a correlation-function anal- 
ysis of a similar dataset. There is a strong motivation for 
performing both analyses. Cosmological measurements from 
large, modern galaxy surveys are limited not by statistical 
errors but by systematic errors, generally arising from mod- 
elling non-linear physics in addition to technical concerns 
such as determining the data covariance matrix in a suf- 
ficiently robust and accurate manner. Systematics imprint 
themselves quite differently into analyses in configuration 
space and Fourier space. For example, a correlation-function 
analysis is more effective at separating processes which oc- 
cur at different spatial scales (for example shot noise ow- 
ing to the discreteness with which galaxies trace dark mat- 
ter haloes is important at all scales in a power-spectrum 
measurement but only at small scales in the corresponding 
correlation function). Moreover, power-spectrum measure- 
ments are much more sensitive to how accurately the se- 
lection function of the survey can be quantified, which is a 
major issue of concern for a survey such as WiggleZ with rel- 
atively low redshift completeness (Blake et al. 2009). How- 



ever, modelling of quasi-linear scales becomes more robust in 
a power-spectrum analysis, where linear theory can be more 
confidently applied across a range of values of wavenumber 
k (whereas structure with a given separation in configura- 
tion space depends on a set of underlying modes with a 
significant range of values of k). In general the consistency 
of the power-spectrum and correlation-function analysis ap- 
proaches is important to check, and would increase confi- 
dence in the robustness of the results against systematic 
error. 

We have investigated in detail the modelling of redshift- 
space distortions in the correlation function in a companion 



paper (iContreras et al. 2011) using dark matter and halo 



catalogues from an N-body simulation. These results are im- 
portant for our current study and we summarize them here: 

• The commonly-used power-law model for the real-space 
correlation function produces a poor fit to the clustering 
pattern and introduces considerable systematic errors in a 
growth rate measurement based on a large sample of galax- 
ies. We considered two alternative models for the real-space 
correlation function, which are able to recover the growth 
rate with minimal systematic error. 

• We introduced a technique which permits a non- 
parametric determination of the galaxy pairwise velocity 
distribution from a 2D correlation function. We used this 
method to measure the pairwise velocity dispersion and 
check the consistency of the assumed functional form. 

• We explored the sensitivity of the results to the range of 
scales included in the fit, in particular considering different 
cuts for the minimum transverse scale fitted, cr m i n . 

Our study of the growth rate in the WiggleZ correlation 
function is structured as follows: in Section 2 we describe our 
data and the measurements of £ s (cr, 7r), and in Section 3 we 
summarize the models fitted. In Section 4 we determine the 
best-fitting model parameters and errors, comparing these 
measurements to the theoretical predictions of the ACDM 
cosmology using standard gravity. We also make an initial 
investigation of the sensitivity of our results to the fiducial 
cosmological model via the Alcock-Paczynski distortion. In 
Section 5 we summarize and discuss the results. 



2 WIGGLEZ SURVEY DATA AND 
CORRELATION FUNCTION 
MEASUREMENTS 

2.1 Data 

The WiggleZ Dark Energy Survey is the latest in a se- 
ries of large spectroscopic galaxy redshift surveys mapping 
out large-scale structure in the Universe. The survey, which 
was carried out between August 2006 and January 2011, 
obtained more than 200,000 galaxy redshifts in the range 
0.1 < z < 1.0 spread across several different equatorial re- 
gions of sky, which together constitute a survey area of 1000 
deg 2 and span a volume of approximately 1 h~ 3 Gpc 3 . The 
observations were carried out at the 3.9m Anglo- Australian 
Telescope (AAT) using the AAOmega spectrograph. In- 
dividual galaxy spectra were redshifted by WiggleZ team 
members. 

A key characteristic of the WiggleZ Survey is the target 



selection criteria ( Blake et al. 2010 Drinkwater et al. 2010) 



which are based on a combination of ultra-violet photome- 
try from the Galaxy Evolution Explorer (GALEX) satellite 
and optical imaging from the SDSS and 2nd Red Cluster Se- 
quence (RCS2) survey ( | Yee et ah) [2005 ) . WiggleZ galaxies 
are star-forming (emission-line) galaxies, which are expected 
to avoid the densest regions of galaxy clusters where star for- 
mation is suppressed by mechanisms such as ram-pressure 
stripping of the gas content of the galaxies. Hence this Wig- 
gleZ dataset should be less susceptible to the non-linear 
clustering effects found on small scales (such as "fingers- 
of-god") compared to surveys targetting more luminous, 
highly-biased galaxies. 

In this study we analyze a total of 162,454 WiggleZ 
galaxy redshifts, which are selected from the [1,3,9,11,15,22]- 
hr regions after minor cuts for redshift and survey contiguity. 
We sub-divide the sample into four redshift slices spanning 
redshifts of [0.1-0.3], [0.3-0.5], [0.5-0.7] and [0.7-0.9]. Table 
[I] shows the details of the catalogues we used. We calcu- 
late an effective redshift z e ff = [0.21,0.39,0.61,0.76] for the 
measurement in each redshift slice by averaging the mean 
redshifts of the galaxy pairs included in our analysis at the 
separations of interest. 



2.2 Random catalogues 

The galaxy correlation function is measured by comparing 
the pair-count of the dataset in separate bins to that of a 
random, unclustered set of galaxies possessing the same se- 
lection function as the data. We generated random WiggleZ 
datasets using the method discussed by |Blake et al.| ( |2010| ). 
This process models several effects including the variation of 
the GALEX target density with dust and exposure time, the 
incompleteness in the current redshift catalogue, the varia- 
tion of that incompleteness imprinted across each 2-degree 
field by constraints on the positioning of fibres and through- 
put variations with fibre position, and the dependence of the 
galaxy redshift distribution on optical magnitude. Random 
catalogues, each containing the same number of galaxies as 
the data, are generated by a Monte Carlo process. 
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^slice/Field 




01 hr 


03 hr 


09 hr 


11 hr 


15 hr 


22 hr 




Zl [0.1-0.3] 


0.21 


3058 


2973 


2671 


4025 


3578 


5782 


22087 


Z2 [0.3-0.5] 


0.39 


2983 


3522 


7438 


9581 


10942 


7728 


42194 


Z3 [0.5-0.7] 


0.61 


5702 


5982 


11294 


13964 


17928 


8924 


63794 


Z4 [0.7-0.9] 


0.76 


3518 


4315 


5189 


7077 


8943 


5337 


34379 


Total 




15261 


16792 


26592 


34647 


41391 


27771 


162454 



Table 1. This table summarizes the number of galaxies used in our analysis, and the effective redshift of our measurement, by region 
and redshift slice. The effective redshift was measured by averaging the mean redshift of each galaxy pair [z m = (z\ + 22)/2] for all the 
galaxy pairs counted in our 2-point correlation function measurements. For the scales of interest (0 — 50 /i — 1 Mpc) we find no significant 
dependence of z e ff on the pair separation. 



2.3 Correlation function measurement 

In order to measure £ s (cr, tt) we first converted the red- 
shifts and celestial coordinates of the galaxies into co-moving 
(x,y, z) spatial co-ordinates assuming a fiducial flat ACDM 
cosmological model with Q m — 0.27 (and Ho — 100 h km 
s _1 Mpc -1 ). We investigate the sensitivity of our results to 



the fiducial model in Section |4.6| We then computed the 
number of galaxy pairs in a series of (a, tt) bins, where tt 
and a are the pair separations parallel and transverse to the 
line-of-sight respectively, where we re-define the line-of-sight 
direction for every pair of galaxies using the bisector of their 
subtended angle on the celestial sphere, in the same man- 



ner as |Hamilton (1993). We repeated this process using 100 
random galaxy catalogues, and finally employed the Landy 
& Szalay estimator for £ s (cr, 7r) (Landy & Szalay 1993): 

DD - 2DR + RR 



RR 



(2) 



where DD is the data-data pair count, DR is the data- 
random pair count, and RR is the random-random pair 
count. We computed these quantities in square (<7,7r) bins 
with side 3/i -1 Mpc over the ranges < a < 30 /i -1 Mpc, 
< tt < 30 h' 1 Mpc. We do not use any additional weights 
for each data or random object, noting that in a sample 
with relatively low number density such as WiggleZ, weight- 
ing results in a negligible improvement to the error in the 
correlation function over the range of scales of interest. The 
fibre-collision corrections are not important for the scales 
analyzed here, and are ameliorated by repeat observations 
of each patch of sky to build up the target density. We mea- 
sured £ s (a, tt) in each separate WiggleZ region and redshift 
slice, and optimally combined the measurements for com- 
mon redshift slices using inverse- variance weighting (see be- 
low). Figure [I] illustrates the combined correlation function 
data for each of the four redshift slices. 



2.4 Covariance Matrix 

We used a jack-knife procedure to measure the correlation 
function variances and covariances between bins. This tech- 
nique is implemented by sub-dividing the dataset into Njk 
equal volumes, then repeating the £ s (a, tt) computation Njk 
times excluding each of the subsets in turn. Using this pro- 
cedure we obtain Njk different, correlated measurements 
of £s(cr, tt). In order to compute the full covariance matrix 
between the bins we use the expression 



(3) 



In our analysis we used Njk = 125 for each individual sur- 
vey region, choosing equal- volume sub-regions tiled in three 



dimensions, corresponding to a total Njk — 750 for the 
combination of the 6 regions. We checked that our results 
were not sensitive to these choices. Excluding the diagonal 
elements, the correlation coefficients from the combined co- 
variance matrices have consistently low average values 0.10 
± 0.10; 0.09 =1= 0.08; 0.07 ± 0.08 and 0.03 ± 0.07 for the four 
redshift slices in ascending redshift order. 

The combination of £ s (cr, tt) for the 6 sky regions for ev- 
ery redshift slice is produced using inverse- variance weight- 
ing: 



(4) 



for the i th bin, and the covariance matrices are combined as 



k=1 a lk a lk 



6 1 6 1 



(5) 



for the i th and the j th bin of £ s (cr, tt). We determined the 
bin size (3/i _1 Mpc) of our measurement as a compromise 
between mapping the (cr, tt) plane with high resolution, and 
the need to obtain an invertible covariance matrix. A smaller 
bin size would require a larger number of jack-knife sub- 
volumes to ensure a non-singular covariance matrix, result- 
ing in smaller jack-knife regions which would be less inde- 
pendent on the scales of interest. 



3 MODELS FOR THE 2D CORRELATION 
FUNCTION 

The effect of linear redshift-space distortions on the power 
spectrum in Fourier space was described by Kaiser (1987): 



Ps(k,n) 



(l+Ptl 2 ) 2 Pr{k) 



(6) 



where P s and P r are the redshift-space and real-space galaxy 
power spectra at wavenumbers k, ji is the cosine of the angle 
of the Fourier mode to the line-of-sight, and ft = f /b quan- 
tifies the amplitude of redshift-space distortions in terms of 
the growth rate / and linear bias factor b. Hamilton (1992) 
provided the equivalent expression in configuration-space: 



tt) = £o(r)Po(M) + &(r)ffeO*) + U(r)P^) 



(7) 



where Pi are the Legendre polynomials and the correlation 
function multipoles £e(r) are given by the general expres- 



The growth rate in WiggleZ 5 



Log 10 (Correlation Function) 




-20 o 

o [Mpc/h] 



Log 10 (Correlation Function) 



i 



-0.5 



- -1,0 



- -1,5 
I 




I 







a [Mpc/h] 




a [Mpc/h] 
Log 10 (Correlation Function) 



a [Mpc/h] 







- -1,0 



2,0 




I 



Figure 1. Measurements of the redshift-space correlation function for redshift slices 0.1 — 0.3 (top left), 0.3 — 0.5 (top right), 0.5 — 0.7 
(bottom left) and 0.7 — 0.9 (bottom right), obtained by combining results in the different WiggleZ Survey regions with inverse- variance 
weighting. Only the top-right quadrant of data for each redshift is independent; the other three quadrants are mirrors of this first 
quadrant. Noticeable is the lack of "fingers-of-god" prominent in similar measurements for LRGs, owing to the avoidance by WiggleZ 
galaxies of high-density regions. 
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In order to complete the model, the 2D correlation function 
including linear redshift-space distortions, £s(<7, 7r), is then 
convolved with a function F(v) representing a dispersion of 



pairwise velocities on small scales: 



f 

J — c 



H(z)a(z) 



F(v) dv 



(13) 



normalized such that F(v) dv — 1. 

The most common approach in previous analyses of 
redshift-space distortions in the galaxy correlation function 
has been to use a power-law choice for £ r and an exponen- 



tial ID velocity distribution for F(v) (Landy & Szalay 1998 



Hawkins et al. |2003 Cabre & Gaztanaga 2009). In a com- 



panion paper (|Contreras et al. 2011) we presented a de- 



tailed analysis of the systematic errors resulting from these 
assumptions, considering three different models for the un- 
derlying real-space correlation function £ r : 

• Power-law model: £ r (r) = 

• CAMB model: we use a real-space correlation function 
calculated numerically for a given set of cosmological param- 
eters (motivated by analyses of the Cosmic Microwave Back- 
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ground radiation) using the CAMB software (Lewis et al. 
2000 ) with the WMAP-7 best-fit cosmological parameters. 

• Quadratic Correlation Function (QCF) model: we gen- 
eralize the power-law function to an empirical real-space 
correlation with one further degree of freedom: £ r (r) = 

J+qlog 10 



This form has more flexibility to de- 
scribe the real-space correlation function than the other 
models. For q = the QCF recovers the power-law, while 
for q ^ —0.5 the QCF is similar to a CAMB correlation 
function. 



We note two advantages of the QCF model. Firstly it 
does not assume a particular set of cosmological parameters, 
so represents a more flexible, model-independent approach 
that is able to reproduce the shape of the CAMB correlation 
function for a range of cosmological parameters. Secondly, 
the flexibility of the QCF model allows it to accommodate 
effects such as scale-dependent bias, which is expected to be 
important for massive haloes hosting Luminous Red Galax- 
ies. 

In Contreras et al. ( 20111) we also considered a series 



of options for the shape of the pairwise velocity distribution 
function F(v), introducing a parameter x which interpolates 
between exponential and Gaussian shapes for this function: 



F(y)=xF(v) e + (l-x)F(y)g 



(14) 



The motivation for adding this parameter was to allow the 
data to determine which of the exponential or Gaussian dis- 
tributions provided the best fit. We found that, when fit- 
ting to halo catalogues from dark matter simulations, the 
small-scale non-linear anisotropy in the clustering pattern 
favoured the exponential form, whilst if we excluded the 
small-scale part of the data by restricting the fitted range to 
o~ > (Jmin the choice of exponential or Gaussian was rendered 
unimportant. In our study of N-body simulations we were 
able to measure the shape of the pairwise velocity distribu- 
tion directly, and the results were self-consistent with the 
best-fitting model for F(v) when the CAMB or QCF mod- 
els were used for the real-space correlation function. When 
the power-law model was used for £ r (r) instead, we found 
systematically higher values for the pairwise velocity disper- 
sion. This result provided more evidence that using an ac- 
curate model for the real-space correlation function yielded 
constraints for all the fitted variables with lower systematic 
error. 

We reformulated the problem of obtaining the pair- 
wise velocity distribution by the process of fitting a non- 
parametric stepwise function instead of a exponential or 
Gaussian profile. We introduced a fast technique to do this 



(Contreras et al. 2011 appendix A). The results for the 



pairwise velocity dispersion measured as a stepwise function 
from the simulation were consistent with both the direct 
measurements, and with the fits of parameterized models. 

Finally, we concluded from our experiments with sim- 
ulations that our models were still not complex enough 
to model the small-scale clustering pattern of dark matter 
halo catalogues, and we were required to introduce a cut 
o~ > cr m in ~ 6/i -1 Mpc in the fitted range to restrict sys- 
tematic errors in the derived growth rate to a maximum of 
Af « 0.05. 



4 RESULTS 



4.1 Fitting procedure 



We fitted the models described in Section [3] to the WiggleZ 
dataset by minimizing the x 2 statistic calculated as 



Y,{M l -D t ) T C- j i {M 3 -D 3 ) 



(15) 



where Mi is the model vector for i data bins, Di is the 
data vector, and is the element of the inverse covari- 
ance matrix which relates bins i and j. Our default redshift- 
space QCF model contained 6 parameters [ro, 7, q, /3, o~ v , x] 
and we considered variations as described in the text. We 
used a Monte Carlo Markov Chain (MCMC) code to per- 
form these fits, determining the best-fitting values and joint 
likelihood distributions of each parameter. We found that 
chains of 100,000 iterations gave robust results for all mod- 
els. We double- checked our results using downhill simplex 



algorithms. Following Contreras et al. (2011), we investi- 



gated the dependence of our results on the minimum trans- 
verse scale a included in our fits. The correlation function 
at low a may contain the signature of strong non-linearities 
not fully described by our models, potentially introducing a 
systematic error in the growth measurement. 



4.2 Measurements of the redshift-space distortion 
parameter f3 

In Figure [2] we illustrate the measurements of f3 in each 
redshift slice, combining results from the different WiggleZ 
regions as a function of the minimum transverse separation 
c"min • We show results for each of the three real-space correla- 
tion function models introduced above. We found signatures 
of systematic errors at small scales with the same pattern 
as found in our investigation of mock halo catalogues from 
N-body simulations, especially when the power-law model 
was assumed. 



4.3 Measurements of the galaxy bias factor b 

We determined the linear galaxy bias factor 6, which should 
be combined with the measurement of f3 to produce the 
growth rate f = /3 b, assuming the CAMB correlation func- 
tion model describes well the dark matter correlation func- 
tion £ m . Then, the bias is determined simultaneously with 
the other parameters when fitting the CAMB model to the 
galaxy correlation function £ g using: 



(16) 



Figure ^illustrates the best-fitting bias in each redshift slice 
as a function of cr m in- 

The results for x and a v show similar behaviour as 
found in our simulation study: for small cr m i n the value of x 
favoured an exponential shape for F(v) and a v was tightly 
constrained, and for larger cr m i n both parameters were not 
well-determined. See Figure [3] for an example. 

The reduced x 2 values of the best-fitting models are 
displayed in Figure [5] We note that the QCF and CAMB 
models both produce good fits to the data, with the QCF 
model carrying the additional advantage of not assuming a 
small-scale shape for the correlation function. 
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Figure 2. Measurements of the redshift-space distortion param- 
eter (3 for the four WiggleZ redshift slices, fitting to the combined 
correlation function of the survey regions. We show results for 
the three real-space correlation function models discussed in the 
text: QCF, CAMB and power-law, as a function of the minimum 
transverse scale cr m i n included in the fit. Note the convergence 
of the three correlation-function models after small-scale bins are 
excluded from the fits. 



The fitted galaxy bias factors must be corrected for the 
(relatively small) effects of redshift blunders caused by the 
mis-identification of emission lines in the WiggleZ survey 
scattering the true position of galaxies. This effect has been 

Section 3.4) and 



carefully quantified by Blake et al. (2010 



causes a loss of signal in the correlation function which may 
be recovered by multiplying by a redshift-dependent correc- 
tion factor: 



: m [1 



(17) 



where r is the fraction of redshift blunders at redshift z. 
The relevant correction factors (1 — r) -2 for the different 
redshift slices can be found in Table [5] and are equivalent to 
an upward correction in the measured galaxy bias factors. 

4.4 Stepwise velocity distribution 

Our stepwise fitting technique allows us to measure the pair- 
wise velocity distribution in a non-parametric fashion, in- 
stead of fixing its shape to some functional form. The 2D 
correlation function of WiggleZ galaxies does not contain 
strong "fingers-of-god" features such as those encountered 
in redder galaxy samples such as SDSS LRGs. This is ex- 
pected, as WiggleZ galaxies are blue star-forming galaxies, 
which inhabit less dense environments than their red coun- 
terparts and are thereby less affected by non-linear clus- 
tering. Figure [6] shows the non-parametric stepwise velocity 
distribution obtained in every redshift slice of our sample. 




2 4 6 
a m in [MpC/h] 



8 



Figure 3. Example of the behaviour of the parameter x as a 
function of the minimum transverse scale cr m i n included in the 
fit, in this case for the third redshift slice. Values of x close to 1 
implies that F(y) has an exponential shape rather than a Gaus- 
sian shape. It becomes less well determined when excluding the 
small-scale part of the data. 



We used seven bins in velocity of width ~ 300 km s~ up 
to a maximum of 2000 km s _1 . We find a slow growth in 
the width of the pairwise velocity distribution as redshift 
decreases, as expected due to the growth of non-linear struc- 
ture with time. Where the data permits us to discriminate, 
the measured stepwise distribution is consistent with an ex- 
ponential function for F(v). 



4.5 Measurements of the growth rate / 

The final growth rate measurements are obtained by com- 
bining f3 and the bias: f(z) = /3 b . However, as the bias 
was measured assuming a CAMB matter correlation func- 
tion with a fixed normalization as , the measured value of b is 
degenerate with as . We can convert our measurements to be 
independent of this effect by quoting the value of f(z)ag(z) 
(Son g &; Percivalj |2009| ). In Figure [7] we display these re- 
sults in the four redshift slices for different choices of the 
fitting range cr m i n . Motivated by the simulation results pre- 



sented by |Contreras et aL (2011), we take as our preferred 
measurement the QCF model with cr m i n = 6/i -1 Mpc. 

In Figure [8] we compare these preferred measurements 
to the results from the analysis of the 2D WiggleZ survey 
power spectra (|Blake et al. 2011a), together with some pre- 



vious galaxy surveys, and overplot the prediction of a ACDM 
model with matter density Q m = 0.27 and two different val- 
ues of the low-redshift normalization ag — 0.7 and 0.8. The 
separate measurements using different datasets and tech- 
niques are consistent within the statistical errors. The er- 
rors in the growth rate measurements using the 2D WiggleZ 
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Figure 4. Measurements of the galaxy bias factor in the four 
WiggleZ redshift slices as a function of the minimum transverse 
scale cr m i n included in the fit. The bias is determined by fitting 
a CAMB model to the 2D correlation function. We note that 
the bias measurements presented in this Figure have not been 
corrected for the effect of redshift blunders (see text). 



Figure 6. The best-fitting stepwise values of the pairwise velocity 
distribution of WiggleZ galaxies. We find that the amplitude of 
the pairwise velocities systematically decreases with increasing 
redshift, which is expected in a model where non-linear structure 
grows with time. 



correlation function are around twice the size of those de- 
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Figure 5. Reduced x 2 values corresponding to the best-fitting 
parameters of each of the correlation function models in the four 
WiggleZ redshift slices, as a function of the minimum transverse 
scale. 



2011a) 



termined by fits to the power spectra (Blake et al. 
We attribute this to a combination of the greater number 
of parameters that must be varied to obtain a good fit to 
the correlation function given the greater relative impor- 
tance of the non-linear effects that must be modelled, and 
the more limited range of scales (cr, tt) < 30/i _1 Mpc stud- 
ied. Nonetheless, the statistical agreement of these analyses 
is evidence that the different measurements of the growth 
rate are not dominated by systematic errors. 



4.6 Effects of the Alcock-Paczynski distortion 

Finally we tested the sensitivity of our results to the fidu- 
cial cosmological model adopted in our analysis by vary- 
ing the value of the matter density fi m . The sensitivity of 
our measurements to the fiducial model is a result of the 



Alcock-Paczynski effect (Alcock & Paczynski 1979), a geo- 



metrical distortion of the measured clustering pattern which 
arises if the trial cosmology differs from the true cosmology. 
The effect is partially degenerate with redshift-space distor- 



tions (Ballinger et al. 


1996 Matsubara & Suto 


1996| Mat- 


subara| |2000| |Seo & ] 


ilisenstein, 2003; Simpson 


h Peacock, 


2010|) resulting in systematic variations in the best-fitting 



growth rate as the fiducial cosmology is altered. Alterna- 
tively, assuming that the cosmological model is known from 
a ACDM standard gravity model with parameters tuned to 
match the Cosmic Microwave Background anisotropics, the 
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Figure 7. The red circles, green triangles and blue diamonds 
represent the fitted values of the growth rate in the four Wig- 
gleZ redshift slices using a QCF, CAMB and power-law model 
respectively for different data cuts a > cr m i n . This plot shows 
the convergence of the three models when more of the small-scale 
part of the data is excluded. We note that our simulation study 
|Contreras et al.| ( |2011| shows that the growth rate is systemati- 
cally under-estimated if cr m i n is too low, with cr m i n ~ 6h~ 1 Mpc 
required to obtain systematic-free results with the QCF model. 
In all panels we overplot the prediction of a fiducial ACDM model 
with Q m = 0.27 and erg = 0.8. 



Alcock-Paczynski effect provides a further cross-check of this 
model. 

We implemented this analysis by repeating our fitting 
procedure for different values of the fiducial matter den- 
sity Q m = (0.07,0.17,0.37,0.47) for the central redshift 
slice z = 0.6. The results are shown in Figure [9] demon- 
strating that the ACDM WMAP-cosmological model is self- 
consistent in the range 0.2 < Q m < 0.3 and produces incon- 
sistent results (in standard gravity) for other values of Q m . 
We note however that our treatment here is approximate; in 
our power spectrum model we have fixed the values of the 
other cosmological parameters at the WMAP7 best-fitting 
values. A full treatment of these model fits is beyond the 
scope of the current study and is provided by Parkinson et 
al. (in prep.). 



5 CONCLUSIONS 

We have analysed redshift-space distortions in the 2-point 
galaxy correlation function of the WiggleZ Dark Energy Sur- 
vey built from over 160,000 galaxies, dividing the dataset 
into four redshift slices. We fitted a series of different mod- 
els to the observed £ s (cr, ty) (2D 2-point correlation function) 





1.0 : 




0.8 : 


s 




"CO 


0.6 \ 






s 






0.4 : 




0.2 ' : 




0.0 - 



• P, Spectrum 
C. Function 



o 8 =0.8 
o 8 =0.7 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 8. Growth rate results from power-spectrum and 
correlation-function analyses (green and red circles) of WiggleZ 
data. In black, we overplot results from the 2dFGRS, SDSS and 
VVDS surveys (in increasing redshift order). Both WiggleZ sets of 
results are shifted in redshift for clarity. We show predictions for 
ACDM cosmological models assuming Q m = 0.27 and as = 0.7 
(dashed line) and 0.8 (solid line). 
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Figure 9. A self-consistency test of the ACDM cosmology as- 
suming standard gravity, in which the WiggleZ growth rate mea- 
surements in the z = 0.6 redshift slice are repeated for different 
fiducial choices of f2 m . The results are sensitive to the Alcock- 
Paczynski effect, and produce a measured growth rate inconsis- 
tent with the prediction of standard gravity for both low and high 
values of Q m . 
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using a Monte Carlo Markov Chain process. We summarize 
our results as follows: 

• We determine the cosmic growth rate across the redshift 
range 0.1 < z < 0.9 with an accuracy of around 20% in each 
redshift slice of width Az — 0.2, and our measurements 
are consistent with the prediction of a ACDM model with 
Q m « 0.27 and a 8 « 0.8. 

• Our results are also consistent with an independent 
power-spectrum analysis of a similar dataset (| Blake et al 



nomical Observatory in the development and support of the 
AAOmega spectrograph, and the running of the A AT. 

This research was supported by C A ASTRO: 



2011b), demonstrating that systematic errors are not signif- 



icant in these two different approaches. 

• We perform a non-parametric determination of the pair- 
wise velocity distribution of WiggleZ galaxies as a step- 
wise function, demonstrating that the amplitude of pairwise 
velocities grows with decreasing redshift, as expected due 
to the growth of non-linear structure with time, and that 
(where it is possible to fit a model) it is well-described by 
an exponential function. 

• Our measurements agree well with the behaviour found 



in our previous experiments using N-body simulations ( Con- 



treras et al. 2011), and we obtain a convergence in the dif- 



ferent models for the growth rate when excluding the small- 
scale non- linear part of the data (a < 6/i _1 Mpc). 

• We repeated the entire procedure (from transform- 
ing redshift catalogues into 3D spatial galaxy distributions 
to running MCMC correlation function model-fitting) for 
ACDM cosmologies with different values of Q m . The results 
show that the ACDM model is only self-consistent in the 
range 0.2 < Q m < 0.3. 

We conclude that the standard ACDM cosmological 
model provides a good description of the growth rate which 
drives redshift-space distortions in the clustering of WiggleZ 
galaxies, using cosmological parameters which also yield a 
good simultaneous fit to the expansion history measured 
by baryon acoustic oscillations and supernovae, and to the 
fluctuations in the Cosmic Microwave Background radiation. 
We have shown that different statistical analyses of the Wig- 
gleZ dataset, with very different potential systematic errors, 
produce growth rate measurements in close mutual agree- 
ment. 
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